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Abstract 

The process e"'"e~ — > 77(7) is studied using data recorded with the OPAL detector at 
LEP. The data sample corresponds to a total integrated luminosity of 25.38 pb^^ taken at 
centre-of-mass energies of 130 - 172 GeV. The measured cross-sections agree well with the 
expectation from QED. In a combined fit using data from all centre-of-mass energies, the 
angular distribution is used to obtain improved limits on the cut-off parameters: A_|_ > 195 
GeV and A_ > 210 GeV (95% CL). In addition, limits on non-standard e"'"e~7 couplings 
and contact interactions, as well as a 95% CL mass limit for an excited electron, M^* > 194 
GeV for an e~'"e~7 coupling k = 1, are determined. 
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1 Introduction 



This paper reports a study of the annihilation process e+e~ 77(7) using data recorded with 
the OPAL detector at LEP. At LEP energies, this is one of the few processes having negligible 
contributions from the weak interaction. Since the QED differential cross-section is precisely 
predicted in theory, deviations from the expected angular distribution are a sensitive test for 
non-standard physics processes contributing to these photonic final states. 

The OPAL collaboration has previously published a study of photonic final states, with and 
without missing energy, at ^/s = 130 - 140 GeV The present analysis concentrates on final 
states with two or more detected photons, but no missing transverse momentum, to study only 
the QED process. Photonic final states with missing energy have been analysed separately 0. 

Any non-QED effects described by the general framework of effective Lagrangian theory should 
increase with centre-of-mass energy. Existing OPAL limits on deviations from QED can be 
improved by using the data at centre-of-mass energies of 161.3 GeV and 172.1 GeV. A small 
amount of data taken at 170.3 GeV is included in the 172 GeV sample. The corresponding 
integrated luminosities of these data sets are 9.97 and 10.13 pb~^, respectively. Since the 
selection criteria have changed, previously analysed data taken at centre-of-mass energies of 
130.3 GeV (2.69 pb^^) and 136.2 GeV (2.59 pb~^) are reanalysed here to allow for a coherent 
treatment. The 136 GeV sample includes a small amount of data taken at 140.2 GeV. The 
error on the luminosity differs slightly for the different energies and is approximately 0.5% . 

These measurements test QED at the highest centre-of-mass energies. Possible deviations are 
conveniently parametrised by cut-off parameters A±. A comparison of the measured photon 
angular distribution with the QED expectation leads to limits on the QED cut-off parameters 
A±, contact interactions (e+e~77) and non-standard e"'"e~7-couplings as described in section 3. 
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The possible effects of an excited electron, e*, which would also change the angular distribution, 
are investigated. In addition, the possible production of a resonance X via e+e~ X7, followed 
by the decay X — 77, is studied in the invariant mass spectrum of photon pairs in three-photon 
final states. 

The following section contains a brief description of the OPAL detector and the Monte Carlo 
simulated event samples. Section 3 describes the QED differential cross-sections for e^e~ — >■ 
77(7), as well as those from several models containing extensions to QED. In sections 4-6 the 
analysis is described in detail. The results are presented in section 7. 



2 The OPAL detector and Monte Carlo samples 

A detailed description of the OPAL detector can be found in 0. OPAL uses a right-handed 
coordinate system in which the z axis is along the electron beam direction and the x axis 
is horizontal. The polar angle, 6, is measured with respect to the z axis and the azimuthal 
angle, 0, with respect to the x axis. For this analysis the most important detector component 
is the electromagnetic calorimeter (ECAL) which is divided into two parts, the barrel and 
the endcaps. The barrel covers the polar angle range of |cos^| < 0.81 and consists of 9440 
lead-glass blocks. The endcaps cover the polar angle range of 0.81 < | cos 6*1 < 0.98 and 
consist of 1132 blocks. In this analysis, the central tracking detector is used primarily to reject 
events inconsistent with purely photonic final states. Raw hit information from the vertex drift 
chamber (CV) and jet drift chamber (CJ) is used to reject events with tracks coming from the 
interaction point. CV is divided into 36 0-sectors and its inner 12 (6) axial layers cover an 
angular range of | cos^| < 0.95(0.97). CJ is divided into 24 ^-sectors and covers an angular 
range of | cos^^| < 0.97 with its inner 16 layers. Incorporated in the surrounding magnet yoke is 
the hadronic calorimeter (HCAL) covering 97% of the solid angle. The outermost detectors are 
the muon chambers, shielded from the interaction point by at least 1.3 m of iron and covering 
the polar angle range of | cos 6*1 < 0.985. 

Different Monte Carlo samples are used to study efficiency and background. For the signal 
process e+e^ — > 77(7) the RADCOR 0] generator is used. It provides a 0{a^) cross section 
up to I cos^l = 1 for the photon angle. No Monte Carlo with complete fourth order is currently 
available. The only program that generates four-photon final states neglects the mass of the 
electron and therefore does not correctly include photons in the far forward range. The Bhabha 
process is studied with two different programs. BHWIDE ^ generates both electron and 
positron in the acceptance of the detector. In contrast, TEEGG allows one of them to have 
very low energy or escape along the beam-pipe, in addition, one photon is scattered into the 
detector. The process e"'"e" Vwy^j) is studied with NUNUGPV 0. Both e"'"e~ /i"*"/!" 
and e"'"e~ —>■ t^t^ are simulated using KORALZ |H| and PYTHIA [3 is used for hadronic 



events. All samples were processed through the OPAL detector simulation program |T^ and 
reconstructed as for real data. 
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3 Cross section for the process e+e — > 77 



The differential cross-section for tlie process e"'"e 
QED is given by [n|: 

I da\ a 

\'dn 



77 in tfie relativistic limit of lowest order 

2 1 + COS2 



Born 



s 1 — cos^ 6 



(1) 



where s denotes the square of the centre-of-mass energy, a is the electromagnetic coupling 
constant and the polar angle of one photon. Since the two photons cannot be distinguished 
the event angle is defined such that cos 9 is positive. 



In Ref. |jT2[ possible deviations from the QED cross-section for Bhabha and M0ller scattering 
are parametrized in terms of cut-off parameters. These parameters correspond to a short range 
exponential term added to the Coulomb potential. This ansatz leads to a modification of the 
photon angular distribution as given in Eq. 

r,2 





1 ± 



2Ai 



sin^^ 



(2) 



Alternatively, in terms of effective Lagrangian theory, a gauge invariant operator may be added 
to QED. Depending on the dimension of the operator different deviations from QED can be 
formulated |13[. Contact interactions (77e"'"e~) or non-standard 7e"'"e~ couplings described by 
dimension 6, 7 or 8 operators lead to angular distributions with different mass scales A (see 
Eqs. 0- 1). The subscripts (QED+6 etc.) follow the notation in Ref. ||T3[ . 





Born 



1 + 
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1 



sin^^ 



.^/Born'^32^Af 



+ 



Born 



s^ml 1 
327r a! 



(3) 
(4) 
(5) 



The definition of Eq. is identical to the standard definition (Eq. ^ if A^ = f Ag. Similarly 
Eq. @) is equivalent to Eq. (^ if Ag = ml A^. Therefore only the parameters of Eq. @) and 
(0) are determined by a fit to obtain limits on deviations from QED. The limits on the other 
parameters can easily be derived from these results. 



The existence of an excited electron e* with an e*e7 coupling would contribute to the photon 
production process via t-channel exchange. The resulting deviation from 

(^)b depends on 



a 



the e* mass Me* and the coupling constant n of the e*e7 vertex [|T^ 
da 



da \ 

.^/Born^ 



+ 4 
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/ K 



Me. 



(e^ sin^ e + Ml 
Ml sin2 6 



Xq^-M^Y {q'^-Ml)\ 



{q^ - M!,){q'^ - MX 
+ 



q^ 



q' 
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q'- 



+ sin^ e 
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+ 
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-.12 



Ml 



with the beam energy E = y/s/2, = —2E'^{1 — cos6) and q'^ = —2E'^{1 + cos0). In the hmit 
Me* ^ a/s, the mass is related to the cut-off parameter by M^* = a/k A+. 



4 Event angle definition and radiative corrections 

For the process e+e" 7172 the polar angle 9 of the event is defined by the angle between 
either of the two photons and the beam direction since | cos ^1 1 = | cos 6*2 1 . This is a good 
approximation for most of the events under consideration, since additional photons tend to be 
soft. For many events, however, there is a third energetic photon and thus | cos6'i| 7^ | cos ^2 1 
in general. Several angle definitions are possible to characterize an event. The following two 
are considered: 

I cos6'i| + I cos6'2| 

COSe'av = 7: , (7) 



cosr 



. 9i — 02 

sm 



sm — - — , (8j 



where 61 and 62 are the polar angles of the most energetic photons. Both cos6'av and cos^* 
are identical to | cos^| for two-photon final states. For three-photon events in which the third 
photon is along the beam direction, 6* is equivalent to the scattering angle in the centre-of-mass 
system of the two observed photons. 

Fig. |1] shows the ratio of the angular distributions using both cos^av and cos^*, relative to 
the Born cross section as derived using an 0{a^) e~^e~ 77(7) Monte Carlo generator 0]. 
The event angles are calculated from the two photons with the highest generated energy. The 
comparison is made at the generator level, i.e. without detector simulation and efficiency effects. 
It can be seen that the distribution of cos6'av (Eq. |^) shows large deviations from the lowest 
order (Born) distribution for much of the cos^ range. For this analysis cos^* (Eq. (§)) is 
chosen because it better matches the shape of the Born distribution over the range cos^* < 0.9 
considered in this analysis. 



5 Event selection 



Events are selected by requiring two or more clusters in the electromagnetic calorimeter (ECAL). 
A cluster is selected as a photon candidate if it is within the polar angle range | cos6'| < 0.97. 
The cluster must consist of at least two lead-glass blocks, with a combined ECAL energy de- 
posit exceeding 1 GeV uncorrected for possible energy loss in the material before the ECAL. 
Events with a photon candidate having five or more reconstructed clusters within a cone with 
a half-angle of 11.5° are rejected. This isolation criterion helps to reduce some instrumental 
background. 

There are two major classes of background remaining to the 77(7) signature. The first can be 
identified by the presence of primary charged tracks. Bhabha events, for example, have similar 
electromagnetic cluster characteristics as 77(7) events, but are normally easily distinguished 
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by the presence of charged tracks. The second class consists of events without primary charged 
tracks. Certain cosmic ray events and the Standard Model process e+e" — > Du^^ contribute to 
this background. 



5.1 Neutral events 

Events having only photons in the final state are classified as 'neutral events'. They should 
not have any charged track consistent with coming from the interaction point. The rejection of 
all events having tracks in the central tracking chambers CV or CJ would lead to an efficiency 
loss because of converted photons. Nevertheless, contributions from any channel with primary 
charged tracks should be reduced to a negligible level. 

To reject events with primary charged tracks while retaining efficiency for converted photons, 
only the inner part of the drift chambers are considered. First, the correlation between the 
observed clusters and charged hit activity in both drift chambers is used. Hits are counted in 
the 0-sectors of CV and CJ which are geometrically associated to each cluster. A correlation is 
assigned to a cluster if there are more than a given number of wires with hits in the associated 
0-sector. 

• A CV correlation is assigned if there are at least m wires with hits in the n CV layers 
nearest to the beam-pipe (denoted by m/n), depending on cos 6* of the cluster: 



Cut on m/n 


cos 9 region 


6/12 


0. 


< 


cos^ 


< 


0.75 


5/8 


0.75 


< 


cos 9 


< 


0.95 


4/6 or 5/8 


0.95 


< 


cos 9 


< 


0.97 



• A CJ correlation is assigned if there are at least 12 wires with hits in the inner 16 C J 
layers, independent of the cluster polar angle. 

Two vetoes are defined using combinations of these hit activity correlations in CV and CJ. A 
third veto tests for reconstructed charged tracks not correlated with either of the clusters. Any 
of the three vetoes rejects the event. 

• The single veto requires that both the CV and CJ correlation are assigned for any 
cluster. 

• The double veto requires that for each of the highest energy clusters either the CV or 
CJ correlation is assigned. 

• The unassociated track veto requires that there be no reconstructed track with a 
transverse momentum of more than 1 GeV and at least 20 hits in CJ, separated by more 
than 10° in from all photon candidates. 



7 



5.2 Cosmic ray events 



A cosmic ray particle can pass through the hadronic and electromagnetic calorimeters without 
necessarily producing a reconstructed track in the central tracking chambers. Events of this 
type are rejected if there are 3 or more hits in the muon chambers. In the case of 1 or 2 muon 
hits the event is rejected if the highest energy HCAL cluster with at least 1 GeV is separated 
from each of the photon candidates by more than 10° in cf). Events are rejected if the cluster 
extent in cos 6 is larger than 0.4. This cut is primarily to reject beam halo events. 



5.3 Kinematic selection 



The event sample is divided into three classes /, // and ///. The classes are distinguished by 
the number of photon candidates and the acoUinearity angle defined as C = 180° — ^, where 
^ is the angle between the two highest energy clusters. Different selections are applied to each 
class separately, to make use of the different kinematics. Only events with cos^* < 0.9 are 
selected to avoid systematic errors due to large efficiency and radiative corrections. 

All events having an acoUinearity angle C, < 10° (i.e. the two highest energy clusters are 
almost coUinear) belong to class / independent of the number of photon candidates. For true 
e+e~ — >■ 77(7) events in this class, the sum of the two highest cluster energies Eg = Ei + E2 
should almost be equal to the centre-of-mass energy ^/s. The distribution of Es is shown in 
Fig. ||. Events having Es > 0.6^/s are selected. This cut is well below the tail of the energy 
distribution for e^e~ 77(7) Monte Carlo events as shown in the figure. 

Class // contains acollinear events {( > 10°) with exactly two observed photon candidates. 
Events of this class typically contain an energetic photon that escapes detection near the beam- 
pipe (I cos 6*1 > 0.97). If the polar angle of this photon is approximated as | cos^^| = 1, its energy, 
£"1081 (Eq. can be estimated from the angles of the observed photons 6*1 and 6*2. The energy 
sum Es is then defined as the sum of the two observed cluster energies and the lost energy. 

^ / sin di + sin 62 \ , . 

\ I sm(6'i + 6'2)|/ 

Es = El+E2 + £lost (10) 

The imbalance B, defined as 



B = (sin^i + sin^2l 



cos 



(11) 



provides an approximate measure of the scaled transverse momentum of the event without 
using the cluster energies. Fig. ^ shows the distributions of B and Es- It can be seen that the 
background (mainly Z/Z777) is uniformly distributed in B whereas the signal is peaked at low 
values. Events are selected if i3 < 0.2 and Es > 0.6y/s. Since the angular definition discussed 
in section H uses the two highest energy photons, events are rejected if i^iost exceeds the energy 
of either observed photon. 

Class /// contains acollinear events > 10°) having 3 or more observed photon candidates. To 
calculate the transverse and longitudinal momenta {pt, pi) of the system the cluster energies 
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have to be used in addition to the photon angles. Since a non-zero longitudinal momentum 
could correspond to an additional photon along the beam direction, the energy sum Es is 
calculated as sum of the cluster energies Ei and pi. 

n 

Es = Y.Er+Pi. (12) 

i=l 

Fig. H shows the distribution of Es/y/s versus Pt/y/s for e"'"e~ — >• 77(7) Monte Carlo and for 
the data. In the data the e+e~ 77(7) events are clearly separated from the background by 
the fact that they have small transverse momenta and an energy sum around the centre-of-mass 
energy. The main part of the background originates from cosmic ray events without hits in the 
muon chambers. The selection requirements Es > O.Q^/s and pt < O.ly^ easily reject these 
events. 

No event with more than three clusters is observed. The angle sum is used to identify 
planar three-photon events: 

^ a = aij + aik + akj, (13) 

with aij the angle between photons i and j. Seven planar events with J2<^ > 350° are accepted 
as three-photon events and included in the sample of e"'"e~ — 77(7). Two other events are 
consistent with three detected photons and an additional photon along the beam direction. 

The kinematic requirements used to select each of the event classes are summarised in Tab. 



Event class 


Requirements 


all 


cosO* < 0.9 


/ 


C < 10° 
Es > 0.6v/^ 


// 


C > 10° 

Es > 0.6^ 
Ei,E2 > -Eiost 

B < 0.2 
2 photon candidates 


III 


C > 10° 

Es > 0.6y^ 

Pt < O.lv^ 
> 3 photon candidates 


planar 


Ea > 350° 


nonplanar 


'Ea < 350° 



Table 1: Summary of the kinematic cuts. For definition of the variables see the text. 



6 Corrections and systematic errors 

Since the deviations from QED (Eqs. |^ - §[) are given with respect to Born level, the observed 
angular distributions need to be corrected to Born level. The effect of radiative corrections to the 
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Born level calculation is quantified by 7^, the ratio of the angular distribution of e+e llil) 
Monte Carlo and the Born cross-section as shown in Fig. 0: 

The ratio TZ is used to correct the data bin by bin to the Born level. A 1% error on the total 
cross-section from higher order effects is assumed. No error on the slope of the distribution is 
included in the results. Since the 0{a^) radiative corrections are small, the 0{a^) effects are 
assumed to be negligible. 

The efficiency and angular resolution of the reconstruction is determined using a Monte Carlo 
sample with full detector simulation. The efficiency is reasonably constant for cos^* < 0.9, but 
drops rapidly for cos^^* > 0.9. The overall efficiency for cos 6** < 0.9 is 91.9% with a maximum 
of 95% in the barrel of the detector. A polynomial parametrisation S{cos6*) is used for the 
efficiency correction. Due to uncertainties of the photon conversion probability and the Monte 
Carlo statistics a 1% systematic error is assumed for the efficiency. The agreement between 
generated and reconstructed angles is very good. An angular resolution of 0.3° full width at 
half maximum is obtained. Background is studied using Monte Carlo events from the processes 
shown in Tab. 0. The expected ratio of background to signal is less than 0.4% and is neglected. 



Process Generator 


Background events 


e+e" ^ e+e" (BHWIDE) 


< 0.24 


e+e- ^ e+e- (TEEGG) 


< 0.67 


e+e" — s> z/z/77 


< 0.02 


e+e" 


< 0.04 


e+e~ — > r+r~ 


< 0.05 


e+e~ — * qq 


< 0.02 



Table 2: Estimated 95 % CL upper limits for expected background processes from Monte Carlo 
at = 130 - 172 GeV. 

The probability that a signal event is rejected by the neutral event selection due to random 
instrumental background causing a veto is studied with randomly-triggered events. For the 
single veto the probability is 4 x 10~^ and it is 1 x 10"'^ for both the double veto and the track 
veto for. The small overall veto probability of 5 x 10~^ is therefore neglected. The systematic 
errors on the total cross-section are summarized in Tab. ^ 





Uncertainty 


Luminosity 

Radiative correction IZ 
Selection efficiency £ 
Background 


0.5% 
1.0% 
1.0% 
< 0.4% 


Total 


1.6% 



Table 3: Summary of systematic errors on the cross-section 
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7 Results 



In Tab. ^ the numbers of observed events in each class are compared to the QED expectations. 
The derived total cross-section a in the range cos^* < 0.9 is plotted in Fig. ^ as a function 
of the centre-of-mass energy. The numerical results for the cross-section are given in Tab. |^. 
They are corrected for efficiency loss and 0{a^) effects. All numbers agree well with QED 
expectations. 



The measured differential cross-sections at 130, 136, 161 and 172 GeV centre-of-mass energies 
are shown in Figs. ^ and |^ together with a fit of the function (^)^ (Eq- 0)- The fit to the 

distribution is performed using the binned log likelihood method. The likelihood function C is 
based on Poisson statistics and defined as: 

rii 

= ^e-^* (15) 
nil 

with Hi the number of observed and Hi the number of expected events per cos 6* bin i. To 
determine /ij the model dependent cross-section function ^ is not integrated over the bin. 
Instead a simple procedure is applied in which the central value Xj of the bin is determined as 
defined in Ref. |T3[ : 

^ ^ {y)dy, (16) 




Xu - Xi J XI \d^lj 



Born 



where xi and Xu are the lower and upper boundaries of the bin i. In this way the differential 
function ^ can be directly compared to the integrated number of events presented as a his- 
togram. The mean efficiency St and radiative corrections TZi are included in the expectation /ij 



Energy [GeV] 


130 


136 




Expected 


Observed 


Expected 


Observed 


Class / 


34.3±0.9 


33 


30.1±0.8 


26 


Class // 


4.0±0.3 


2 


3.5±0.3 


3 


Class /// planar 


1.0±0.2 


2 


0.9±0.2 





Class /// nonplanar 












^Born 
"tot 


15.7 


14.9±2.5 


14.3 


12.2±2.2 


Energy y/s [GeV] 


161 


172 




Expected 


Observed 


Expected 


Observed 


Class / 


85.8±0.6 


90 


77.6±0.7 


75 


Class // 


9.0±0.2 


8 


8.0±0.2 


14 


Class /// planar 


1.9±0.1 


3 


1.8±0.1 


2 


Class in nonplanar 




1 




1 


^Born 
"tot 


10.2 


10.9±1.1 


9.0 


9.7±1.0 



Table 4: Comparison of number of observed events and Monte Carlo prediction. For nonplanar 
events no expectation is given, since the 0{a^) Monte Carlo does not include these events. 
The two observed class /// nonplanar events are kinematically compatible with a fourth photon 
along the beam-direction. In addition the total cross section corrected to the Born level is 
given. 
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for each bin. To allow the total number of expected events to vary within the systematic error, 
a normalization factor e is added: 

yUi = e ^(^i) i^u - Xi)£iniL (17) 

1 rii 

s^ = - J2s{cose^), (18) 

where L is the integrated luminosity and cos^j is the angle of the j-th event. An estimator 
function P is defined which includes a Gaussian term with mean 1 and width 5 = 0.016 (see 



Tab. 0) to account for the error of the normalization e. The routine MINOS |16], which provides 
asymmetric errors, is used to minimize P: 



52 

-I? 



52 +E2(Aii- n^ln/ii). (19) 



The fit is performed with two free parameters: the normalization e and the model dependent 
parameter A (see Tab. ^ and Eqs. ^ | and |^). To obtain the limits at 95% confidence level the 
probability is normalized to the physically allowed region, i.e. A+ > and A_ < as described 
in Ref. IT^ 



Results for the different parameters are obtained from a simultaneous fit to the angular disti- 
butions for each centre-of-mass energy. A fit is also performed for each centre-of-mass energy 
separately, with the results for A± given as an example in Tab. ^. The limits for the combined 
fit are summarised in Tab. To determine the limit on the mass of an excited electron Mg* a 
fit is performed using f^j ^ (Eq. ^. For the results given in Tab. |^ the couphng constant k 



Model 



do_\ 



A 



s [GeV] 



±1/A1 



± 



130 
136 
161 
172 
130 - 172 



(da\ 



\l/A.^\ 130-172 



\dn)e* 



130 - 172 



Fit result 



A 



/4 3+20.6\ 
V "^-18.0^ 
y 1 +19.3 
' -16.7 

+5.12 
-4.70 



1.53 



10-10 GeV-^ 
10-10 GeV-^ 
10-10 GeV~4 
10-10 GeV-^ 
10-10 GeV-4 



5.57+J;J) ■ 10-18 GeV-6 



8.4+11;°) . GeV-2 



0.999 ±0.016 
0.999 ±0.016 
1.001 ±0.016 
1.000 ±0.016 
1.000 ±0.016 



1.000 ±0.016 



1.000 ±0.016 



Table 5: Results for fit parameters A and e. For A± the results for all energies are shown 
seperately. The error on the normalisation e refects the assumed systematic error. 
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for the (e*e7)-vertex is fixed at k = 1. Fig. ^ shows the upper hmit (95 % CL) on versus 
the mass of an excited electron Mg*. 



The angular distributions for all energies agree well with the QED expectation. The lower 
limits obtained from the combined fit on A±, A and M^* are higher than existing published 
results using lower energies (see Ref. |T^, [EDI and PT|) and are in agreement with other 



results obtained at this centre-of-mass energy [22|. Previous limits on the excited electron mass 
with K = 1 are Me* > 129 GeV f^TJ, 136 GeV and 147 GeV 0. 



A resonance X produced in the process e+e" X7 and decaying photonically X ^ 77 would 
be seen in the two photon invariant mass spectrum since this process leads to a three-photon 
final state without missing energy. This search has been performed previously at the Z'^-peak 
23 1 and at higher energies |21]. The invariant mass of each photon pair is shown in Fig. ^ for 
all events of classes // and ///. There are three entries for events with three clusters. Since 
the angular resolution is very precise, the energies of the three photons are calculated from the 
angles assuming three photon kinematics: 



Ek oc sin aij ; Ei + E2 + E-^ = a/s, 



(20) 



with Ek the energy of one photon and a^j the angle between the two other photons. For 
class // events |cos^| = 1 is assumed for the unobserved photon. A typical mass resolution 
for photon pairs of about 0.5 (0.7) GeV can be achieved for class /// (//). The distribution 
agrees well with the Monte Carlo expectation from the QED process e"'"e~ —>■ 77(7), with no 
enhancement due to a resonance is observed. From the class /// distribution an upper limit 
on the total production cross-section times the photonic branching ratio of an isotropically 



produced resonance is calculated using the method of Bock ||2^. Combining the data of all 
centre-of-mass energies and subtracting the e^e~ — > 77(7) background the limits shown in Fig. 
p!0| are obtained. The mass range is defined by the phase space of the selection and limited due 
to the acoUinearity restriction. 



Parameter 


A+ 


A_ 


Ae 


At 


As 


Me* 


[GeV] 


195 


210 


793 


483 


15.5 


194 



Table 6: Summary of 95% CL lower limits obtained from the combined fit to the ^/s = 130, 
136, 161 and 172 GeV angular distributions. The results are for the cut-off parameters A± and 
mass scales A according to QED+6, QED+7 and QED+8 expectation (Eqs. |]- Ag and Ag 
are derived from A+ and A7 respectively. The lower limit for the mass of an excited electron is 
also determined with the coupling constant k, assumed to be = 1. 



8 Conclusions 



The QED process e"'"e~ 77(7) has been studied using data taken with the OPAL detector 
at LEP energies above the Z° resonance. Both the angular distributions and the total cross- 
section measurement agree well with QED predictions. Limits are set on cut-off parameters, 
mass scales for contact interactions (77e"'"e~) and for non-standard 7e"'"e~ couplings, as well 
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as on the mass of an excited electron coupling to e7. These limits are listed in Tab. ^ In 
the 77 invariant mass spectrum of events with three final state photons, no evidence is found 
for a resonance X decaying to 77. No photonic event with four or more detected photons is 
observed. 
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Figure 1: Ratio of the differential cross-section for the e"'"e — > 77(7) Monte Carlo sample 



relative to the Born cross-section, TZ 



Born 



TZ is shown here for both angular 



definitions cos6'av (Eq- 0) and cos6'* (Eq. |). The event angle is calculated from the two 
photons with the highest generated energy. 
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Figure 2: Scaled sum of the two highest cluster energies for all events with an acoUinearity 
angle C < 10° (corresponding to class /). The points with error bars represent the data, the 
histogram the e+e" — > 77(7) Monte Carlo expectation. The cut on this quantity is indicated. 
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Figure 3: Event distributions for data in class II for ^/s = 130, 136, 161 and 172 GeV. Plot 
a) shows the distribution of the imbalance B, a measure of the scaled transverse momentum 
(for definition see Eq. 11) together with the selection cut. Plot b) shows the scaled sum of 
both cluster energies plus Eiost after the cut on B. The cut is indicated. The points represent 
the data, the solid histogram the Monte Carlo expectation from e"'"e~ — > 77(7) and the dashed 
histogram the Monte Carlos expectation from background (mainly Z/I777). 
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Figure 4: The scaled energy sum versus the scaled transverse momentum for class /// events, 
for a) e+e~ — > 77(7) Monte Carlo and b) the OPAL data. The box indicates the selected 
region. The background comes mainly from cosmic ray events. 
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Figure 5: Total cross-section for the process e"'"e~ — > 77 with cos6 < 0.9. The data are corrected 
for efficiency loss and higher order effects and correspond to a Born level measurement. The 
result at the Z° is taken from Ref. ||19|. The curve corresponds to the Born level QED prediction. 
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Figure 6: The measured angular distribution for the process e"'"e~ 77(7) as selected in the 
three classes at ^/s = 130 and 136 GeV. The data points show the efficiency-corrected number 
of events; radiative corrections are also included. The solid curve corresponds to the Born level 
QED prediction. The dotted lines represent 95% CL intervals of the fit to the function ( ^ ) 
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Figure 7: The measured angular distribution for the process e"'"e~ 77(7) as selected in the 
three classes at ^/s = 161 and 172 GeV. The data points show the efficiency-corrected number 
of events; radiative corrections are also included. The solid curve corresponds to the Born level 
QED prediction. The dotted lines represent 95% CL intervals of the fit to the function ( ^ ) 
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Figure 8: Upper limit (95 % CL) on the square of the couphng constant as a function of the 
mass of an excited electron Mp. . 
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Figure 9: The invariant mass of photon pairs for a) class // events and b) class /// events. The 
points are the data, the histogram the e'^e" — > 77(7) Monte Carlo expectation. There is one 
entry per event for class // events and three entries per event for class /// events. 
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Figure 10: Lower limits (95 % CL) for the cross section times branching ratio for the process 
e'*"e~ — > X7, X ^ 77 as a function of the mass of the resonance X. The e"'"e~ 77(7) 
background is subtracted. The step at 120 GeV comes from the phase space hmit of the 130 
and 136 GeV data. 
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